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Abstract. Three-quarters of the 1 km'^ neutrino 
telescope IceCube is currently taking data. Current 
models predict high-energy neutrino emission from 
transient objects Uke supernovae (SNe) and gamma- 
ray bursts (GRBs). To increase the sensitivity to such 
transient objects we have set up an optical follow- 
up program that triggers optical observations on 
multiplets of high-energy muon-neutrinos. We define 
multiplets as a minimum of two muon-neutrinos from 
the same direction (within 4°) that arrive within a 
100 s time window. When this happens, an alert 
is issued to the four ROTSE-III telescopes, which 
immediately observe the corresponding region in 
the sky. Image subtraction is appUed to the optical 
data to find transient objects. In addition, neutrino 
multiplets are investigated online for temporal and 
directional coincidence with gamma-ray satellite ob- 
servations issued over the Gamma-Ray Burst Coor- 
dinate Network. An overview of the full program is 
given, from the online selection of neutrino events to 
the automated follow-up, and the resulting sensitivity 
to transient neutrino sources is presented for the first 
time. 

Keywords: Neutrinos, Supernovae, Gamma-Ray 
Bursts 

I. Introduction 

When completed, the in-ice component of IceCube 
will consist of 4800 digital optical modules (DOMs) 
arranged on 80 strings frozen into the ice, at depths 
ranging from 1450m to 2450m [1]. Furthermore there 
will be six additional strings densely spaced at the 
bottom half of the detector The total instrumented 
volume of IceCube will be 1 km''. Each DOM contains 
a photomultiplier tube and supporting hardware inside a 
glass pressure sphere. The DOMs indirectly detect neu- 
trinos by measuring the Cherenkov light from secondary 
charged particles produced in neutrino-nucleon interac- 
tions. IceCube is most sensitive to neutrinos within an 
energy range of TeV to PeV and is able to reconstruct 
the direction of muon-neutrinos with a precision of ~1°. 
The search for neutrinos of astrophysical origin is among 
the primary goals of the IceCube neutrino telescope. 
Source candidates include galactic objects like super- 
nova remnants as well as extragalactic objects like 
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Fig. 1. Neutrino event spectrum in the IceCube detector, from kaon 
and pion decay in the supernovae-jet model of Ando and Beacom [5]. 

Active Galactic Nuclei and Gamma-Ray Bursts [9] [10]. 
Offline searches for neutrinos in coincidence with GRBs 
have been performed on AMANDA and IceCube data. 
They did not lead to a detection yet, but set upper 
limits to the predicted neutrino flux [13]. While the 
rate of GRBs with ultra-relativistic jets is small, a much 
larger fraction of SNe not associated with GRBs could 
contain mildly relativistic jets. Such mildly relativistic 
jets would become stalled in the outer layers of the 
progenitor star, leading to essentially full absorption of 
the electromagnetic radiation emitted by the jet. Hence, 
with the postulated presence of mildly relativistic jets 
one is confronted with a plausible but difficult-to-test hy- 
pothesis. Neutrinos may reveal the connection between 
GRBs, SNe and relativistic jets. As was recently shown, 
mildly relativistic jets plowing through a star would be 
highly efficient in producing high-energy neutrinos [5]- 
[7]. The predicted neutrino spectrum follows a broken 
power law and Fig.[T]shows the expected signal spectrum 
for neutrinos produced in kaon and pion decay in the 
source, simulated using the full IceCube simulation 
chain. The expected number of signal events is small 
and requires efficient search algorithms to reduce the 
background of atmospheric neutrinos (see section HH). 
An optical follow-up program has been started which 
enhances the sensitivity for detecting high-energy neutri- 
nos from transient sources such as SNe. In this program, 
the direction of neutrinos are reconstructed online, and 
if their multiplicity pass a certain threshold, a Target- 
of-Opportunity (ToO) notice is sent to the ROTSE-III 
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network of robotic telescopes. These telescopes monitor 
the corresponding part of the sky in the subsequent 
hours and days and identify possible transient objects, 
e.g. through detection of rising supernova light-curves 
lasting several days. If in this process a supernova is 
detected optically, one can extrapolate the lightcurve or 
afterglow to obtain the explosion time [2]. For SNe, 
a gain in sensitivity of about a factor of 2-3 can be 
achieved through optical follow-up observations of neu- 
trino multiplets [4] . In addition to the gain in sensitivity, 
the follow-up program offers a chance to identify its 
transient source, be it a SN, GRB or any other transient 
phenomenon. 

II. Neutrino Alert System 

IceCube's optical follow-up program has been operat- 
ing since fall of 2008. In order to match the requirements 
given by limited observing time at the optical telescopes, 
the neutrino candidate selection has been optimized to 
obtain less than about 25 background multiplets per year 
The trigger rate of the 40 string IceCube detector is about 
1000 Hz. The muon filter stream reduces the rate of 
down-going muons created in cosmic ray showers dra- 
matically by hmiting the search region to the Northern 
hemisphere and a narrow belt around the horizon. The 
resulting event stream of 25 Hz is still dominated by 
misreconstructed down-going muons. Selection criteria 
based on on track quality parameters, such as number 
of direct hitfl track length and likelihood of the recon- 
struction, yield a reduced event rate of 1 event/(10 min). 
The optimized selection criteria are relaxed to improve 
the signal efficiency, 50% of the surviving events are 
still misreconstructed down-going muons, while 50% 
are atmospheric neutrinos. During the antarctic summer 
2008/2009, 19 additional strings were deployed, which 
have been included in the data taking since end of 
April 2009. To take into account an enhancement in the 
rate due to the increased detector volume, the selection 
criteria have been adjusted and will yield a cleaner event 
sample containing only 30% misreconstructed muons. 
From this improved event sample, neutrino multiplet 
candidates with a time difference of less than 100 s 
and with an angular difference (or 'space angle') of less 
than 4° are selected. The choice of the time window 
size is motivated by jet penetration times. Gamma-ray 
emission observed from GRBs has a typical length of 
40 s, which roughly corresponds to the duration of 
a highly relativistic jet to penetrate the stellar enve- 
lope. The angular difference is determined by IceCube's 
angular resolution. Assuming single events from the 
same true direction, 75% of all doublets are confined 
to a space angle of 4° after reconstruction. Once a 
multiplet is found, a combined direction is calculated 
as a weighted average of the individual reconstructed 
event directions, with weights derived from the estimated 

'Hits that are measured within [-15ns, 75ns] from the predicted 
arrival time of Cherenkov photons, without scattering, given by the 
track geometry. 



direction resolution of each track. The resolution of the 
combined direction is up to a factor of 1/V2 better 
than that of individual tracks. The multiplet direction 
is sent via the network of Iridium satellites from the 
South Pole to the North, where it gets forwarded to the 
optical telescopes. At this point in time, due to limited 
parallelization of the data processing at the South Pole a 
delay of 8 hours is accumulated. In the near future, the 
online processing pipeline will be upgraded, reducing 
the latency drastically to the order of minutes. 
A total of 14 alerts have passed the selection criteria and 
were sent to the telescopes within 7 months of operation. 

III. Optical Follow-Up Observations 

At the moment IceCube alerts get forwarded to 
the Robotic Optical Transient Search Experiment 
(ROTSE) [3]. Additions to the list of participating tele- 
scopes are planned. ROTSE-III is dedicated to observa- 
tion and detection of optical transients on time scales of 
seconds to days. The original emphasis was on GRBs 
while it more recently has also started a very successful 
SN program. The four ROTSE-III telescopes are in- 
stalled around the world (in Australia, Namibia, the USA 
and Turkey). The ROTSE-III equipment is modest by 
the standards of modern optical astronomy, but the wide 
field of view and the fast response permit measurements 
inaccessible to more conventional instruments. The four 
0.45 m robotic reflecting telescopes are managed by a 
fully-automated system. They have a wide field of view 
(FOV) of 1.85° X 1.85° imaged onto a 2048 x 2048 
CCD, and operate without filters. The cameras have a 
fast readout cycle of 6 s. The limiting magnitude for a 
typical 60 s exposure is around 18.5 mag, which is well- 
suited for a study of GRB afterglows during the first 
hour or longer. The typical full width at half maximum 
(FWHM) of the stellar images is smaller than 2.5 pixels 
(8.1 arcseconds). Note that ROTSE-IIIs FOV matches 
the size of the point spread function of IceCube well. 
Once an IceCube alert is received by one of the tele- 
scopes, the corresponding region of the night sky will 
be observed within seconds. A predefined observation 
program is started: The prompt observation includes 
thirty exposures of 60 seconds length 0. Follow-up 
observations are performed for 14 nights. Eight images 
with 60 seconds exposure time are taken per night. The 
prompt observation is adjusted to the typical rapidly de- 
caying lightcurve of a GRB afterglow, while the follow- 
up observation of 14 days permits the identification of 
an increasing SN lightcurve. Once the images are taken, 
they are automatically processed at the telescope site. 
Once the data is copied from the telescopes, a second 
analysis is performed off-line, combining the images 
from all sites. Image subtraction is performed according 
the methods presented in [8]. Here the images of the 

-Once the delay caused by data processing at the South Pole (see 
section |ll) is reduced to the order of minutes, the prompt observation 
will include ten short observations of 5 seconds, ten observations of 
20 seconds and twenty long exposures of 60 seconds. 
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first night serve as reference, while the images from the 
following nights are used to search for the brightening 
of a SN lightcurve. 

IV. Sensitivity 

The sensitivity of the optical-follow up program is 
determined by both IceCube's sensitivity to high energy 
neutrino multiplets and ROTSE-III's sensitivity to SNe. 
We will distinguish two cases: The first being that no 
optical counterpart is observed over the course of the 
program (assuming 25 alerts per year) and the second 
that a SN is identified in coincidence. 

A. No Optical Counterpart Discovered 

With no coincident SN observed, one obtains an upper 
limit on the average number of SNe that could produce 
a coincidence: iV^Qrpgg < 2.44 (for 90% confidence 
level). Constraints on a given model are obtained by de- 
manding that the model does not predict a number in ex- 
cess of the SN event upper limit. We construct a simple 
model based on Ando & Beacom type SNe [5]. We in- 
troduce two parameters: The first is the rate of SNe pro- 
ducing neutrinos p — (4/37r)^^10^''p2e^_4Mpc^'^yr^^. 
Note that ^26^-4 — 1 corresponds to one SN per year 
in a 10 Mpc sphere, about the rate of all core-collapse 
SNe in the local Universe [12]. Since we expect only 
a subset of SNe to produce high energy emission, one 
can assume P2e^_4 < 1. The second parameter is the 
hadronic jet energy Ej^t = 3 • lO^^eggg^ergs and we 
choose to scale the flux normalization of the model 
of Ando & Beacom, Fq, by £32*51. Fig. |2] shows the 
constraints that one can place on the density and jet 
kinetic energy in the Ejct — P plane. The basic shape of 
the constraints that can be obtained in the E-^ct — P plane 
can be understood from the following considerations. 
The number of neutrinos depends on the jet energy and 
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duce two neutrinos. Therefore the number of detection 

^ic/ROTSE is given by iVfc%oTSE oc P2?-4 
For normalization we use Ando & Beacom-like SNe, 
which occur at a rate of ^26^-4 — 1 with GRB-like 
energies (£3651 = 1) and yield ^f(^^']^oxSE ~ "^^^ 
expected IceCube/ROTSE coincidences per year. 
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According to [11] a non-detection limits the number of 
IceCube/ROTSE coincidences at a 90% confidence level 
to ^ic^ROTSE ^ 2.44. Using Eq. [T] one obtains the 
two-dimensional constraints on density and hadronic jet 
energy for this model: 
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Fig. 2. Sensitivity in the Bjot —p plane after one year of operation of 
the 40 string IceCube detector (dashed line — 90% CL; solid line — one 
coincident detection per year). 



which is a reasonably good representation of the two 
dimensional constraints for not too small densities 
P2?-4 Z 10"^- For GRB-like energies (e^^r^^ = 1), it 
follows that at most one out of 80 SNe produces Ando & 
Beacom-like jets in its core. Phrased in absolute terms, if 
no SN will be detected, the rate of SNe with a mildly rel- 
ativistic jet should not exceed p = 3.1 • 10~^Mpc^'^yr~^ 
(at 90% confidence level) in our program. The cut-off at 
small densities visible in Fig. |2] is due to ROTSE-III's 
limiting magnitude. The sphere (i.e. effective volume) 
within which ROTSE-III can detect Supernovae has a 
radius of about 200-300 Mpc. ROTSE-III effectively 
cannot probe SN subclasses that occur less then once 
per year within this sphere. 

B. Significance in case of a detection 

Next we address the case that a SN was detected 
in the follow-up observations. The task mainly consists 
of computing the significance of the coincidence. We 
compute this for one year of data and 25 alerts. Each 
alert leads to the observation of a Afl = 1.85° x 1.85° — 
3.4 square degree field, hence over the course of the 
year ROTSE-III covers a fraction of the sky given 
by An/An x A'^aiorts = 2.1 • lO^^^ ^g^t assume that 
the time window for a coincident of an optical SN 
detection and candidate neutrino multipet is given by 
Atd, the accuracy with which we can determine the 
initial time of the supernova explosion. Studying the 
lightcurve of supernova SN2008D, which has a known 
explosion start-time given by an initial x-ray flash, we 
have developed an accurate way to estimate Atd from 
a SN lightcurve [2]. We fit the light curve data to a 
model that postulates a phase of blackbody emission 
followed by a phase dominated by pure expansion of the 
luminous shell. Explosion times can be determined from 
the lightcurve with an accuracy of less than 4 hours. A 
detailed description of this method can be found in [2]. 
The number of accidental SNe found will be propor- 
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tional to Atd and the total number of SNe per year that 
ROTSE-III would have sensitivity to detect, if surveying 
the sky at all times, iVnoTSE ~ 10"*. Putting all this 
together the number of random coincidences is: 
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An Atd Atd 

A^alcrtsA^ROTSE^ X — = 0.056—^. (3) 

47r yr a 



For A'^bg ^ 1 this corresponds to the chance probability 
p — 1 — exp(— A^bg) ~ -^bg of observing at least one 
random background event. For Atd = Id and no other 
information, the observation of a SN in coincidence with 
a neutrino signal would have a significance of about 2a. 
The significance can be improved by adding neutrino 
timing information as well as the distance information 
of the object found. We first discuss the extra timing 
information. So far we have only required that two 
neutrinos arrive within 100 s to produce an alert. Thus, 
in the analysis presented above, the significance for two 
events 1 s apart would be the same as for 99 s difference. 
Since the probability pt to find a time difference less 
than At^ due to a background fluctuation is given by 
Pt — Aty/lOOs assuming a uniform background, we in- 
clude the time difference in the chance probability. Next 
we discuss the use of the SN distance. One can safely 
assume that there will be a strong preference for nearer 
SNe, since these are most likely to lead to a neutrino 
flux large enough to produce a multiplet in IceCube. 
Using the distance dsN as an additional parameter one 
can compute the probability to observe a background SN 
at a distance d < dsN- The probability is given by the 
ratio of SNe observed by ROTSE-III within the sphere 
dsN to ah SNe: pd = A'^ROTSE(rf)/^ROTSE- In case of 
a detection both dsN and At^, will be available. We use 
a simple Monte Carlo to obtain the significance of this 
detection. For example the detection of two neutrinos 
with a temporal difference of =10 s in coincidence 
with a SN in dsN =20 Mpc distance has a p-value of 
5 ■ 10~^, which corresponds to S.Sti, assuming a total of 
-^alerts = 25 alcrts found in the period of one year. 

V. Coincidences with GCN-GRBs 

According to current models, about every 15-20th 
GRB that can be detected by IceCube will produce 
a neutrino doublet. Hence there is a small possibility 
that we will find a doublet in coincidence with a GCN 
alert, a case that we consider separately here. The 
significance of such a coincidence can be estimated with 
calculation analogous to Eq. |3] The number of accidental 
coincidences with a time difference less than At is given 
by: 

An At „Ai 

iValcrtsiVcCN^ X — = 3.2 • IQ-^ — . (4) 

47r yr Is. 
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where we have assumed 200 GCN notices and 30 
multiplets a year. A coincidence occurs whenever the 
neutrinos and the GRB overlap within predefined win- 
dows in direction and time. For illustrative purposes, if 
we choose a 1.5-degree directional window and a 4- 
hour time window (corresponding roughly to IceCube's 



point spread function and to GRB observations and 
modeling), Eq. |4] yields an expected background count 
of Nbg = 4.7 • 10^^. This corresponds to a 3.5cr 
effect, or equivalently the expectation of a false positive 
from background once every 2100 years. We can further 
reduce the expected background by assuming that the 
neutrino signal is most likely to be emitted at the same 
time as the gamma rays. Since the background multiplets 
will be distributed uniformly across the 4-hour window, 
we can multiply the chance probability above by the 
factor 

tORB — ti. 
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where the absolute value is taken since we assume 
the neutrinos are equally likely to be emitted before 
the gamma-rays as they are after Note that our flat 
probability assumption for the relative emission times of 
gamma rays and neutrinos from GRBs can, of course, be 
modified to follow any particular theoretical model. With 
all these assumptions, if we observe a coincidence that 
is 300 seconds from the GRB onset time, the chance 
probability is then given by Nbg ■ Pt — 4.7 • 10^** • 
300/14400 = 9.8 • 10"^ which corresponds to a 4.4cr 
result. 

VI. Conclusion 

We have presented the setup and performance of 
IceCube's optical follow-up program, which was started 
in October 2008. The program increases IceCube's sen- 
sitivity to transient sources such as SNe and GRBs and 
furthermore allows the immediate identification of the 
source. Non-detection of an optical counterpart allows 
the calculation of a limit on model parameters such as 
jet energy and density of SN accompanied by jets. 
In addition multiplets of neutrinos are tested for coinci- 
dences with GCN messages. Even a single coincidence 
detection would be significant. 
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